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ABSTRACT 

We present a new approach in modelling the polarized Galactic synchrotron emission 
in the microwave range (20-100 GHz), where this radiation is expected to play the 
leading role in contaminating the Cosmic Microwave Background (CMB) data. Our 
method is based on real surveys and aims at providing the real spatial distributions of 
both polarized intensity and polarization angles. Its main features are the modelling 
of a polarization horizon to determine the polarized intensity and the use of starlight 
optical data to model the polarization angle pattern. Our results are consistent with 
several existing data, and our template is virtually free from Faraday rotation effects 
as required at frequencies in the cosmological window. 

Key words: polarization. Galaxy, cosmic microwave background. Method: numerical. 



1 INTRODUCTION 

■ The polarized component of the diffuse background emis- 
sion in the microwave range is of great interest for both the 
Galactic structure and the CMB. Actually, its measurement 

■ leads to probing the structure of the interstellar medium 
(ISM) and the Galactic magnetic field. Moreover, the detec- 

. tion of CMB Polarization (CMBP) allows the investigation 

' of the early Universe. 

CMB anisotropies and polarization are powerful tools 
to determine cosmological parameters (Sazhin & Benitez 
1995, Zaldarriaga, Spergel & Seljak 1997, Kamionkowski & 
Kosowsky 1998). However, although anisotropies have been 
already detected and space missions (MAP*, PLANCKt) 
are expected to make all-sky surveys down to 0°.l angu- 
lar resolution, CMBP still represents a challenge for as- 
tronomers. The first detection has been just claimed by 
DASI (Kovac et al. 2002) and several experiments will ad- 
dress it soon (SPOrt'f (see Carretti et al. 2002, Cortiglioni et 
al. 2002), MAP, PLANCK, B2K2 (Masi et al., 2002), BaR- 
SPOrt (Zannoni et al., 2002) and AMiBA (Kesteven et al., 
2002) among the others). 

Besides the CMBP low emission level (3-4 on sub- 
degree scales and < 1 on large ones) , difficulties in its 

* http://map.gsfc.nasa.gov/ 

t http:/ /astro. estec.esa.nl/SA-general/Projects/Planck/ 
t http://sport.bo.iasf.cnr.it 



detection are mainly related to the presence of foreground 
noise from Galactic and extragalactic sources. Extragalac- 
tic foregrounds essentially consist of radio and infrared dis- 
crete sources, whereas Galactic foregrounds are generated by 
synchrotron, free-free, thermal dust and spinning/magnetic 
dust emissions. 

Synchrotron polarized emission should represent the 
most relevant foreground in the microwave range: free-free is 
fainter (< 4 at 30 GHz in total intensity, see Reynolds & 
Haffner 2000) and almost unpolarized, whereas thermal dust 
has a polarization degree much smaller than synchrotron 
(Prunet et al. 1998, Tegmark et al. 2000). Evidence for spin- 
ning or magnetic dust emission has been found (Kogut et al. 
1996, de Oliveira-Costa et al., 2002) but it seems to play an 
important role only up to ~ 50 GHz. Moreover, it should 
have a low polarization degree (Lazarian & Prunet 2002). 

In spite of its importance, synchrotron emission is 
scarcely surveyed: existing data mainly cover the Galactic 
Plane area at frequencies up to 2.7 GHz, far away from the 
cosmological window (Duncan et al. 1997, hereafter D97, 
Duncan et al. 1999, hereafter D99, Uyaniker et al. 1999, 
Gaensler et al. 2001, Landecker et al. 2002). The Leiden 
data (Brouw & Spoelstra 1976, hereafter BS76) cover high 
Galactic latitudes, but are limited to < 1.4 GHz and are 
largely undersampled. This situation makes having a reliable 
synchrotron polarized emission template in the 20-100 GHz 
range very important. This would allow, for instance, re- 
liable numerical simulations to set-up and test destriping 
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techniques or foreground separation methods (Revenu et al. 
2000, Sbarra et al. 2003, Tegmark ct al. 2000 and references 
therein). At present, only toy models exist, which do not 
account for the real spatial distribution of both polarization 
intensity and polarization angles (Kogut & Hinshaw 2000, 
Giardino et al. 2002). 

In this paper we present a new approach in modelling 
the Galactic diffuse synchrotron polarized emission in the 
20-100 GHz range. It is based on real surveys and fitted to 
the real spatial distribution of both polarized intensity and 
polarization angles. Low frequency data are used to model 
the polarized intensity and optical starlight is used to model 
polarization angles. This allows the construction of Q and U 
maps covering about half of the sky with the SPOrt experi- 
ment angular resolution (FWHM = 7°). Although our work 
is finalized to the SPOrt experiment, the method is general 
enough to be suitable also for smaller angular scales as soon 
as complete sets of data with sub-degree angular resolution 
will be available. 

The great advantage of this new approach is to produce 
Q and U maps free from Faraday rotation, allowing a direct 
extrapolation to the cosmological window. 

The outline of the paper is as follows: the synchrotron 
polarized emission model and the procedure to build the 
template are presented in Section 2, results and compar- 
isons with existing data are described in Section 3 and 4, 
respectively, whereas Section 5 contains the conclusions. 



2 THE MODEL 
2.1 Ingredients 

The aim of our work is to generate template maps of the 
two linear Stokes parameters Q and U of the Galactic syn- 
chrotron polarized radiation in the cosmological window 
near 100 GHz with the SPOrt angular resolution (FWHM 
= 7°). We divide the problem in two parts: 

(i) constructing a polarized intensity (Jp) map: it can be 
obtained from existing total intensity (7) sky surveys as- 
suming a model linking the polarized to the total intensity 
synchrotron emission; 

(ii) building a map of polarization angles not affected 
by the effects of Faraday rotation. At present, only optical 
starlight data fulfil this requirement. 

The Haslam map (Haslam et al. 1982) is the most complete 
sky survey at radio wavelenghts, where synchrotron emission 
is dominant. It is a full-sky map at 408 MHz with a reso- 
lution of 51 arcmin obtained combining observations taken 
with different radiotclescopes. However, it is not perfect for 
our aims since the free- free emission is still significant, expe- 
cially in the Galactic plane (Reich & Reich 1988, hereafter 
RR88). Consequently, identification and subtraction of this 
contribution is mandatory. As described in Section 2.2, we 
perform this separation using the Dodelson (1997) formal- 
ism, which requires a second map at different frequency. We 
use the Reich (1982) map at 1.4 GHz with an angular resolu- 
tion of 35 arcmin, the only other available survey with abso- 
lute calibration covering a large part of the sky {5 > —17°). 

The polarization angles are taken from the Heiles 
starlight polarization catalogue (Heiles 2000). These data 



have many advantages with respect to those in the radio 
band: they are free from Faraday rotation and they cover 
almost all the sky. The catalogue lists polarization degree, 
position angle and distance of about 9000 stars from both 
hemispheres. Further properties of this catalogue and con- 
siderations confirming its validity as a polarization angle 
pattern for our model are described in section 2.4. 



2.2 Synchrotron Intensity Map 

The different frequency behaviour of free-free and syn- 
chrotron emissions allows the application of the Dodelson 
technique (Dodelson 1997) to Haslam (0.408 GHz) and Re- 
ich (1.4 GHz) maps. The original Dodelson formalism is 
centred on the CMB frequency dependence, so that here 
a slight modification is introduced to adjust the method to 
the synchrotron-free-free case. For each pixel i a vector Ti is 
defined whose elements are the pixel antenna temperatures 
at the two frequencies. Its expression in terms of synchrotron 
and free-free components is given by: 



T^H-Tf-fN,, (1) 

where Tf and Tj^ are the synchrotron and free-free contri- 
butions, respectively, and Ni is the noise. 

Assuming the noise from the two maps is completely 
uncorrelated (it comes from different experiments), the cor- 
relation matrix 



Cab,i = (Na,iNb,i) 



(2) 



is diagonal with pixel variances as elements (o and b indicate 
the frequencies). 

Provided the frequency behaviours (shapes) Ff and F^'-'^ 
are known, the estimators 0| and of Tf and T^-'^ can 
be expressed as: 

©I = 0tFt 

(Aff — fjff-pff 

i i i 

Clff . 



(3) 
(4) 



where 61 and are the unknown amplitudes of synchrotron 
and free-free, respectively. 

In the range 0.408 — 1.4 GHz free- free and synchrotron 
are known to follow power laws, so that their shapes are 



= I/"' 
fP. = V-' 



(5) 
(6) 



with the free- free spectral index a = 2.1 (RR88). The fre- 
quency behaviour of synchrotron radiation depends on the 
energy distribution of relativistic electrons and is spatially 
varying across the sky. 

We define a scalar product between vectors as 



(7) 



This expression is similar to that of Dodelson, apart from 
the normalization factor 



v^2 jpffpffn-^' 



(8) 



which here is based on the free-free shape rather than on 
that of CMB. Again following Dodelson, the best estimate 
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Figure 1. Synchrotron spectral index map in tlie 0.408-1.4 GHz 
range (see text and RR88 for details). The map is in Galactic 
coordinates with the Galactic Centre in the middle and Galactic 
longitudes increasing toward left. 

of the amplitudes is given by 

et = ^i^^/F^T„ (9) 

= J2^o'Fl-T. with j = s,ff, (10) 

where the matrix K is defined as 

K,j=F'-F'. (11) 

As a result, the two components are separated providing the 
two maps of synchrotron {Of) and free- free (6'/^). 

The spectral index /3i has been modelled using the anal- 
ysis of RR88 who compared data at 408 MHz and 1.4 GHz 
obtaining the following results: 

• /3 ~ 2.85 toward the Galactic anticentre at & = 0° of 
Galactic latitude, with a flattening with incresing z {z is the 
height above the Galactic plane); 

• /3 ~ 3.1 in the inner disk (Galactocentric distance < 
8 kpc and < 2 kpc). In this region the spectral index is 
nearly constant, beginning to decrease from |6| ~ 10° until 
reaching the value of /3 ~ 2.65 at |6| ~ 30°; 

• /3 ~ 2.65 for |6i > 30° independently of the Galactic 
longitude I. 

The transition from the inner disk region to the rest of the 
Galactic plane occurs between I = 45° and / — 55°, where 
a flattening from P — 3.1 to (3 = 2.85 is observed. From 
spectral index profiles (cfr. Figure 5 in RR88) we find the 
linear behaviour to be a good approximation for (3 in the 
transition regions. From these considerations, we model the 
distribution of synchrotron spectral indeces as follows (see 
Figure 1): 

1) a region towards the Galactic centre (|6| < 10° and 
/ < 45°, Z > 315°) with /? = 3.1; 

2) a region towards the Galactic anticentre on the Galactic 
plane (6 = 0° and 55° < Z < 305°) with /? = 2.85; 

3) a region at high Galactic latitude (|&| > 30°) with f3 — 
2.65. 

4) the spectral index follows a linear behaviour in the tran- 
sition regions and continuity is imposed at borders. Figure 2 
shows the spectral index behaviour in three special cases. 
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Figure 2. Spatial behaviour of the synchrotron spectral index /3 
for three special cases: 1 = 0° (a), Z = 180° (b) and 6 = 0° (c). 



2.3 From Total to Polarized Intensity Map 

A common result of radio-surveys in polarization is the iden- 
tification of two main components: a strong emission from 
discrete sources (Supernova Remnants - SNRs - and several 
sources with no /-counterpart) and a weaker, diffuse emis- 
sion from a background component (D97, D99, Landecker 
et al. 2002, Gaensler et al. 2001) that appears to be rather 
constant with the longitude independently of both angu- 
lar resolution and frequency. This isotropic background sug- 
gests the presence of a polarization horizon, a local screen 
beyond which the polarized emission is cancelled out (D97, 
Gaensler at al. 2001, Landecker et al. 2002). The horizon 
can be imagined as a sort of bubble centred in the observer 
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Figure 3. Contributions to the Galactic polarized radiation are 
limited by a polarization horizon, a local bubble of radius ijp^ 
centred on the Sun (see text for details). The thickness L{l,b) 
of the synchrotron emitting region of our simple spherical model 
(the distance of the halo point P from the Sun position) is also 
shown. 



position (see Figure 3): the net polarized signal is only that 
integrated along the line of sight out to the horizon, whereas 
signals beyond the horizon are depolarized by variations of 
polarization angles (changes and turbulence in the Galactic 
magnetic field). 

A further element suggesting the existence of such a 
horizon comes from Landecker et al. (2002), who show that 
only the closest SNRs are well visible also in polarized emis- 
sion, whereas the most distant ones completely disappear. 

The size of the horizon is not yet known: it depends on 
several effects along the line of sight, like Galactic magnetic 
field turbulence and electron density variations. However, it 
has been suggested that it can range from 2 kpc (Gaensler 
et al. 2001) up to 7 kpc (D97, Landecker et al. 2002), so that 
a few kpc appeares to be a quite acceptable estimate. 

The polarization horizon allows us to model the relation 
between polarized and total intensity synchrotron emissions. 
Given the mean total synchrotron emissivity J''{v,l,b) at 
Galactic coordinates {l,b) and the thickness L{l,b) of the 
synchrotron emitting region (see Figure 3) in the same di- 
rection, the brightness temperature T^[v,l,b) at frequency 



r{u,l,b)L{l,b), 



(12) 



where K is the Boltzmann constant. The thickness L{l,b) 
depends on the geometrical model describing the space dis- 
tribution of the relativistic-electron gas responsible for syn- 
chrotron emission. As a first step in modelling the polarized 
synchrotron radiation we consider the simplest case where 
the gas is uniformly distributed in the Galactic halo. This is 
represented by a sphere of radius ii = 15 kpc centred into 
the Galactic Centre (GG). Thus, in our simple case the line 
of sight L{1, b) is the distance between the Sun and the edge 
of this sphere: 



L{l,b) — d cos(6) cos(i) 



1 + 



(i?Vrf2 - 1) 



cos2(6) cos2(i) 



(13) 



where d = 9 kpc is the Sun distance from GC (See Figure 3 
and Appendix A for details). 

The polarized brightness temperature Tp can be simi- 
larly defined provided the emission is integrated out to the 
polarization horizon Rph and a polarization degree p is in- 
troduced: 



T;{u,i,b) = 



pJ"{i^, I, b) Rph- 



(14) 



Finally, equations (12) and (14) provide the relation between 
polarized and total intesity emissions: 



T;{u,i,b)=p 



R 



ph 



L{l,b) 



T'{u,l,b) 



(15) 



The quantity p Rph is unknown and represents a free param- 
eter to be calibrated with real data. 



2.4 Polarization angle map 

The propagation of an electromagnetic wave of wavelength 
A through a plasma in presence of a magnetic field B is 
affected by Faraday rotation. The net effect is a change in 
the polarization angle (j) by 



A(j> = RM 



RM = 812 / rie (cm"^)B(AtG) • dl(kpc)rad m^ 



(16) 



where RM is the rotation measure, Ue is the plasma electron 
density and dl is the infinitesimal path along the line of sight. 

Estimates of RM from extragalactic radio sources give 
typical values ranging from tens to hundreds rad/m^ at 
medium and high Galactic latitudes, and from tens to thou- 
sands rad/m^ in the Galactic plane (Simard-Normandin & 
Kronberg 1980, Sofue & Fujimoto 1983, hereafter SF83, 
Brown & Taylor 2001, hereafter BTOl). In particular, the 
behaviour along the Galactic plane is well fitted by (BTOl) 



RM(;) = RMo cos(/ - lo) 
RMo = -183 ± 14 rad/m^ 
lo = 84° ±4°. 



(17) 



Brown & Taylor suggest that the modulation in RM occurs 
because of a local constant magnetic field. 

Equation (17) gives, in the frequency range of the cos- 
mological window (20 — 100 GHz), negligible Faraday rota- 
tion effects (±2.5° at 20 GHz): all we need is a template of 
intrinsic polarization angles. When used at 2.7 GHz, which is 
the highest frequency of present polarization surveys, equa- 
tion (17) results in angular rotations up to ±160°. This 
means that radio polarization data cannot be used to build 
a reliable template of intrinsic polarization angles. 

To overcome this problem we use the Heiles catalogue 
on starlight polarization, the optical frequency being unaf- 
fected by Faraday rotation. 

The polarization vector of starlight is parallel to the 
Galactic magnetic field B because of selective absorption 
by interstellar dust grains, whose minor axis is aligned with 
B (Fosalba et al., 2001). Since the synchrotron polarization 
vector is perpendicular to B, starlight polarization angles 
can be used as a template provided a 90° rotation is per- 
formed. 

Most of the Heiles catalogue stars (~ 87%) are within 
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(degrees): eacli star is represented by a pixel in a HEALPix format 
map with Nside = 128 (pixel size ~ 0°.5)). The map is in Galactic 
Figure 4. Number of Heiles stars versus Sun distance. coordinates. 



2 kpc (see Figure 4), tracking the local magnetic field. Their 
distance is in the order of the polarization horizon size (see 
Section 2.3) confirming that the Heiles catalogue can be 
safely used as a template for the polarization angles of syn- 
chrotron emission. The main problems to face with when 
using Heiles polarization angles are the irregular distribu- 
tion of data and the variable sampling distance. However, 
if the uniformity scale of the angles is compatible with the 
sampling distance of the catalogue, the lack of data can be 
filled by linear interpolation. 

An estimate of the uniformity scale can be obtained 
from Figures 9 and 10 of D97, showing that background 
emission regions have a polarization angle pattern varying 
slowly on scales of 5° — 10°. Only areas with strong sources 
show a more complex structure, but their modellization is 
out of the purposes of this paper. 

The sampling distance of the Heiles catalogue (Figure 6) 
is compatible with the uniformity scale everywhere but in 
the region centred in {I = 135°, b = 40°) where it is greater 
than 10°. We exclude this region from the interpolation pro- 
cedure as well as from final template maps. 

We perform the interpolation by generating Q, U pairs 
corresponding to the Heiles polarization angles 6: 

Qg = cos(2e); (18) 
Ug = sin (261). (19) 

Then, for each pixel of the template map under construc- 
tion we linearly interpolate the Qg and Ug values of the 
three closest stars and compute the corresponding polariza- 
tion angle. The interpolation methods uses parallel trans- 
port as described in Bruscoli et al. (2002). 

2.5 The Procedure 

The polarized synchrotron emission template is built as fol- 
lows: 

(i) the CMB emission and the absolute calibration error 
are removed from both the Reich and the Haslam maps using 
values suggested in RR88 (3.7 K and 2.8 K for Haslam and 
Reich maps, respectively); 



Sampling distances 



















/ / 1 










wwww 




— i — 




























1 


1 




\ 














1 














1 






















u 

























U.O H 10.(1 ili'i' 



Figure 6. Map of sampling distances (degrees) of the Heiles cat- 
alogue. The map is in Galactic coordinates. 



(ii) the resulting maps are resampled in HEALPix^ for- 
mat with Nside = 128, corresponding to a pixel size of about 
half a degree; 

(iii) Reich map data are smoothed to the Haslam resolu- 
tion (FWHM = 51 arcmin); 

(iv) the technique for component separation described in 
Section 2.2 is applied by using the synchrotron spectral in- 
dex pattern previously described. This results in the I syn- 
chrotron emission map; 

(v) the relation between 7 and Ip described by equa- 
tion (15), which introduces the effects of the polarization 
horizon, is used to provide the shape of the Ip map, the pa- 
rameter pRph being still free. Its calibration is performed 
a posteriori on BS76 data and a detailed discussion is pre- 
sented in Section 3; 

(vi) from the Ip map and from the angle map obtained 
from starlight data, Q and U maps are computed; 

(vii) the Q and U maps produced in this way are con- 
volved with a FWHM — 7° Gaussian filter to match the 
SPOrt angular resolution. The smoothing procedure applies 

§ http:/ /www. eso.org/science/healpix/ 
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the parallel transport method described in Bruscoli et al. 
(2002) 



3 THE POLARIZED SYNCHROTRON 
TEMPLATE 

Our results are shown in Figure 7, where the Q &z U tem- 
plates at 1.4 GHz are presented. Figure 8 shows a compari- 
son between the polarized intensity Ip of our model and the 
Ip map obtained from BS76 data. The comparison is per- 
formed only for Ip because polarization angles are strongly 
affected by Faraday rotation at 1.4 GHz (see Section 2.4). 
Areas not covered by our template arc in grey: they corre- 
spond to the Southern sky not surveyed by Reich and the 
region around the North Celestial Pole where starlight angle 
data are too sparse. 

The well known feature of real data, i.e. that Galactic 
plane and high Galactic latitudes having comparable emis- 
sions, is well reproduced in our template. 

Our template is also able to reproduce the brightest 
structures in BS76 data, namely: 

(i) the Fan region (the region situated in the Galactic 
plane at 120° < I < 150°); 

(ii) the North Galactic Spur; 

though the latter is fainter than in real data. On the other 
hand, in our template a feature appears in the Galactic 
Plane towards I = 30°-40°, which is not present in BS76. 
One reason might be the sparse sampling of BS76 in this 
area, where the sampling distance is ~ 4° versus an angu- 
lar resolution of 0.6°, and a source might well have been 
missed. Another possibility could be due to Faraday depo- 
larization effects. A qualitative analysis of the whole BS76 
data set (0.408-1.4 GHz) suggests that at 1.4 GHz only the 
Fan region and the North Galactic Spur are free from depo- 
larization effects. This seems to be confirmed by Junkes et 
al. (1990) who find the polarized intensity decreasing from 
I ~ 50° towards the Galactic Centre with a relevant min- 
imum around / ~ 30°. They argue this behaviour might 
be due to depolarization effects: in particular, at I ~ 30°, 
thermal material in the foreground (Scutum arm) might be 
responsible for the observed low polarization. 

The good agreement between the map obtained from 
BS76 data and our model allows us to calibrate our tem- 
plate (the parameter p Rph is still free) by matching the Ip 
emission of the two maps in a well defined area. We use 
the Fan region, the most defined and morphologically simi- 
lar area in both the two maps. The calibration is performed 
with the 820 MHz BS76 data (see Figure 2 in Bruscoli et al. 
2002) rather than with those at 1.4 GHz, because of their 
better sampling, providing 

p Rpf^ = 0.9 ± 0.09 kpc, (20) 

where the error is dominated by the ^ 10% uncertainty on 
BS76 data calibration. Assuming the polarization p on 7° is 
in the range ~ 0.15-0.3 (Tegmark et al. 2000) we obtain for 
the polarization horizon: 

3 kpc < J?ph < 6 kpc, (21) 
in good agreement with present estimates (D97, Gaensler 



Table 1. Peak emission {Fan region) and Prms of our template 
at the four SPOrt frequencies. The Prms is computed on the low 
emission areas (the faintest 50% pixels). 



1/ (GHz) Ip peak (/^K) Prms(MK) 



1.4 


5 X IQS 


6.6 X 10* 


22 


130 


17 


32 


43 


5.6 


60 


6.5 


0.84 


90 


1.9 


0.25 



et al. 2001, Landecker et al. 2002). This is a further con- 
firmation provided by our model of the relation between 
synchrotron / and Ip including a polarization horizon. 

We extrapolate the Q and U templates at 1.4 GHz to 
the cosmological window, and in particular to the SPOrt 
frequencies: 22, 32, 60, 90 GHz. We use a power law with the 
mean synchrotron spectral index /3 = 3.0 found by Platania 
et al. (1997) in the 1 — 19 GHz range. We do not show 
the resulting maps being the same at 1.4 GHz apart from 
the normalization. Instead, we report the emission of the 
most important structure (the Fan region) and the mean 
polarization level 

Prms = V{Q^) + {U'')) (22) 

of the low emission areas (the faintest 50% pixels) in Table 1 
for all the SPOrt frequencies. 

4 COMPARISONS WITH EXISTING DATA 
4.1 Power spectra 

As a first check we compute the Angular Power Spectra 
(APS) of our model and compare them with APS obtained 
from BS76 maps at 1.4 GHz (Bruscoli et al. 2002). 

To account for the irregular sky coverage of our tem- 
plate wo use the method described by Sbarra et al. (2003) 
and based on Q and U two-point correlation functions. 

These are estimated directly on our maps as 

(6) = Af Af X = Q,U, Ip (23) 

where Af is the pixel i content of map X, i and j identify 
pixel pairs at distance 6. 

The polarized power spectra Cf'^ and C^^ are obtained 
by integration 

/■TT 

Cf = Wt [C"2(6>)Fi,e(^)-FC'^(6>)F2,«(6')]sin(^)d6' 
Jo 

/•TT 

Cf = We [C'^(6l)Fi,«(e)-FC'2(e)F2,e(e)]sin(0)d6> 
Jo 

Cl" = 2ttWi I C'''{e)Pe{cos0)d9, (24) 
Jo 

where the functions Fi,<m and F2,em are described by Zal- 
darriaga (1998), Pe are the Legendre polynomials, and the 
function 

We = e'^'+'^''" (25) 



© 0000 RAS, MNRAS 000, 1-12 



A New approach for a Galactic Synchrotron Polarized Emission Template in the Microwave Range 7 









































w 


\\ 














\ 


\ \ \ \ \ ' . 




-mi 










1 






1 1 
















1 






1 /' 
















/ 




/// 
















1 /////// 
























-5. 


A 




i 






01 K 

































/// / 












A 


\ \ \ 


• •///// 


/ / ( 


1 

1 








V 














( 1 








r 






















^ — 
J — 


1 






/' 
























////// 




















/ 










H 













Figure 7. Q (left) and V (right) maps of our synchrotron polarized emission template (Kelvin). The maps are at 1.4 GHz convolved 
with a FWHM= 7° Gaussian filter. 



Brouw & Spoelstrci 




Figure 8. Same of Figure 7 but for 7p. Our template (left) is compared with 1.4 GHz BS76 data (right). 



with G — FWHM/\/81n2 accounts for beam smearing ef- 



fects. Finally, the total polarization spectrum 
defined as 



C 



cf + cf. 



is simply 



(26) 



The resulting power spectra show significant fiuctuations 
also at high multipoles (see Figures 9 and 10). We find that 
the errors cr(C^) are significantly smaller than these fluctu- 
ations, suggesting that they are intrinsic. Nevertheless, the 
overall behaviour of the power spectra can be represented 
by power laws 



cT = AyI 



Y = E,B, P, /, 



(27) 

provide the values ay 



Linear fits to the quantities lnC|' 
listed in Table 2. 

Bruscoli et al. (2002) find consistent values in their anal- 
ysis of large portions of the BS76 maps, namely 1.2 < ay < 
2 {la C.L.) in the 10 < ^ < 70 range. 

The angular behaviour of real polarized synchrotron 
emission is thus well reproduced by our template. 



Table 2. Best fit APS slopes and amplitudes as obtained from 
our template at 1.4 GHz in the 2 < £ <20 range. 



ay 



Ay (K^) 



2.1 ± 0.3 
1.6 ±0.2 
1.9 ±0.2 
1.9 ±0.3 



0.020 ±0.01 
0.007 ±0.005 
0.210 ±0.04 

0.008 ± 0.005 



4.2 Free-free Emission Map 

To understand if we are subtracting the right free- free contri- 
bution from low frequency total intensity data, we compare 
our free-free map with the Galactic H II region catalogue of 
Kuchar & Clarke (1997). 

This is an all-sky flux compilation at 4.85 GHz of 
760 objects, representing the most comprehensive H II re- 
gion catalogue to date. However, a quantitative comparison 
is not straightforward because the catalogue does not take 
into account the diffuse component: only overall patterns 
can be compared. In our free- free map (see Figure 11) the 
emission is concentrated on the Galactic plane, in particu- 
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Figure 9. (top) and (bottom) power spectra computed 
from our synchrotron polarized emission template. Best fit curves 
are also shown. 



lar towards both the Galactic Centre and the area between 
75° < / < 90°. Furthermore, it is very low at mid and high 
Galactic latitudes. Figure 11 also shows the distribution of 
the Kuchar & Clarke H II regions: they too are concentrated 
on the Galactic plane and their number is larger where our 
map shows the largest and strongest structures. This is also 
suggested by Figure 12, which compares the free- free emis- 
sion from our map with the H II region number-density along 
the Galactic plane. Our map is well traced by the free-free 
emitting sources, making us confident that we are subtract- 
ing thermal contributions properly. 

A further test is allowed by the study of thermal emis- 
sion along the Galactic plane described in RR88. 

A relevant contribution (up to 15% at 408 MHz) is 
found in the 15° < I < 50° area as well as a concentration of 
H II regions at I ~ 75° in the Cygnus, where a spectral-index 
flattening reveals the presence of a strong thermal compo- 
nent. These areas correspond to the most evident structures 
in our map, showing the template is also able to reproduce 
important free-free diffuse emission sources. 

Finally, a comparison of our free- free map at 1.4 GHz 
with the Reich map at the same frequency shows that in 
some areas on the Galactic plane the estimated free-free 
contribution is 50% or more of the total emission. This is 




10 



Figure 10. The same as Figure 9 but for C'p (top) and 
(bottom) power spectra. 



an a posteriori confirmation that the free-free contribution 
in the Haslam and Reich maps is not negligible. 

4.3 The Impact of Faraday Rotation 

A check of our polarization angle map is performed by com- 
paring Heiles angles with the data of the Parkes Southern 
Galactic Plane survey at 2.4 GHz (D97). 

However, as pointed out in Section 2.4, Parkes data can- 
not be representative of intrinsic polarization angles because 
of their high RM values. To account for them we introduce 
a compensation by using equation (17). As shown by SF83, 
this behaviour is correct in the Parkes area but in the region 
around I ~ 300°, where a strong deviation is observed and 
a typical value of RM cannot be defined. 

The Parkes survey covers the Galactic plane in the I < 
6°, 238° < I < 360° and |6| < 5° area with a resolution of 
~ 10' at 2.4 GHz. 

We smooth both maps on 4° angular scale to limit the 
impact of local RM vaxiations. Larger scales cannot be ad- 
dressed since they are marginally compatible with the 10° 
width of the Parkes survey. Furthermore, as already dis- 
cussed, starlight data are rotated by 90° to match the syn- 
chrotron polarization angles orientation. 

Figure 5 shows that Heiles angles vary smoothly with 
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Figure 11. Map of free-free emission at 1.4 GHz resulting from 
our procedure (top). Map of Galactic H II region number density 
taken from the Kuchar & Clarke catalogue (bottom). 



the Galactic longitude. Parkes data are slowly varying as 
well but in correspondence of extended sources (D97). Here 
data show peculiar features (discontinuities, inversions, sud- 
den rotations) which even the smoothing procedure is not 
able to remove, the sources extending over several degrees. 
These regions are excluded from our comparison since the 
RM model synthesized by equation (17) just describes the 
general behaviour of background emission. In details, fol- 
lowing the D97 identification we exclude the regions 260° < 
/ < 272° (Vela SNR), 272° < / < 285° (a bright source with 
no total intensity counterpart), 320° < / < 340° (SNR), 
0° < Z < 6° (Galactic centre with several peculiar struc- 
tures). 

We divide the rest of the Parkes surveys in six patches, 
of at least 8° x 4° , characterized by a small variation of the 
polarization- angle pattern. 

For each selected patch we average the difference be- 
tween Parkes and Heiles polarization angles, the latter being 
rotated by 90°. Should Heiles angles describe the magnetic 
field responsible for synchrotron emission, these differences 
would match the polarization angle variations induced by 
RM. These two quantities are reported in Figure 13, whereas 
their difference, expected to be zero, is shown in Figure 14. 

The general agreement between Parkes-Heiles differ- 
ences and RM effects confirms that Heiles data provide a 
reliable template for the polarization angles of Galactic syn- 
chrotron emission. 

Finally, we stress that the evident disagreement at ap- 




Galactic longitude I ( ° ) 



Figure 12. Free-free emission at 1.4 GHz along the Galactic 
plane (6 = 0°) as obtained from our procedure (solid). The data 
arc smoothed on 2° angular scale. The number density of the H II 
region of the Kuchar & Clarke catalogue is shown for comparison 
(dashed). 




360 340 320 300 280 260 240 



Figure 13. Differences between polarization angles in the Parkes 
survey and in our template (filled). Faraday rotation values at 
2.4 GHz as resulting from equation (17) are shown for comparison 
(open). 

proximatively I ~ 300° is not surprising, being RM data in 
this area not well fitted by equation (17). 

Both SF83 and D97 note that this region presents a 
complex situation where a sudden inversion of polarization 
angles takes place, probably due to the transition from the 
Carena to the Centaurus arms (/ ~ 302°). Here a 180° ro- 
tation of magnetic field occurs generating large changes in 
RMs. 



5 CONCLUSIONS 

In this paper we have presented a new approach for a tem- 
plate of the polarized Galactic synchrotron emission which, 
free from Faraday rotation effects, can be better extrap- 
olated to the cosmological-window frequency range (20- 
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Figure 14. Differences between the two quantities plotted in 
Figure 13. The dashed line shows the expected value. 

100 GHz). Differing from previous spatial models (Giardino 

et al. 2002, Kogut & Hinshaw 2000), it is intended to pro- 
vide the real spatial distribution of both polarized inten- 
sity and polarizationn angles. We notice that most previous 
works adopted a complementary approach based on angular 
frequency rather than real space (Tucci et al. 2000, 2002; 
Baccigalupi et al. 2001; Giardino ct al. 2001; Bruscoli et al. 
2002). In fact, angular spectra are commonly used for scale 
separation in the case of CMB, since they are suitable for 
cosmological parameters fitting. However, the shape of po- 
larization angular spectra found at frequencies ~ 1 GHz, 
being affected by Faraday rotation, cannot be confidently 
extrapolated to the cosmological window. This point was 
raised up by Tucci et al. (2001) and Bruscoli et al. (2002), 
who noticed different behaviours in Cf and C^'' spectra at 
all scales in the range / ~ 10 -i- lO'*, suggesting the latter be 
less affected by Faraday rotation. In particular, the analy- 
sis of the ATCA Test Region at 1.4 GHz (Tucci et al. 2002) 
shows strong changes of slope at small angular scales for 
with X — E. B. P, but not for Cj^ , and this seems to bo 
the most dramatic effect of Faraday screens. At present, no 
method is known for correcting such effects directly on angu- 
lar spectra. The present model is intended to overcome such 
a problem too: polarization angular spectra in the cosmo- 
logical window should bo computed on the spatial template 
rather than simply extrapolated from the direct analysis of 
low-frequency maps. 

The model construction consists in three steps: 

1) The synchrotron total intensity / is estimated from the 
low-frequency radio surveys, cleaned from free-free emission 
(see Section 2.1 and 2.2). 

2) The polarized intensity Ip is estimated from I by as- 
suming that the Galactic polarized synchrotron emission is 
local because of the existence of a polarization horizon (see 
Section 2.3). 

3) The polarization angle map is built up from starlight 
polarization data (sec Section 2.4). 

Step 3 is of great importance for building a pattern of 
Stokes parameters. The available RM measurements suggest 
in fact that the effects of Faraday rotation on polarization 



angles arc still too relevant in the radio-surveys at 2.7 GHz, 
the highest available frequency, so that the intrinsic posi- 
tion angles cannot be estimated even with RM corrections 
in large portions of the sky. In our approach we simply over- 
come the problem using the starlight optical data (Hciles 
2000). The local origin of this catalogue (~87% of the stars 
within 2 kpc) and its frequency unaffected by Faraday ro- 
tation effects make it a reliable template for polarization 
angle of the Galactic synchrotron. Our analysis shows also 
that the sampling of the catalogue is compatible with the 
SPOrt angular resolution in all the sky but in the North 
Galactic Pole where it is too sparse. 

A set of checks provides the consistency of the model 
with existing data: 

• The free-free map obtained with our procedure well 

traces the HII region distribution from the Kuchar & Clarke 
(1997) catalogue. This makes us confident about the validity 
of step 1. 

• The estimate provided for the distance of the polariza- 
tion horizon (3-6 kpc) is in good agreement with the values 
obtained by observations, and the polarized intensity Ip well 
reproduces the main structures observed in the BS76 data 
at 1.4 GHz. Both facts support the reliability of steps 1 and 
2. 

• The slopes of polarized angular power spectra Cf , Cf , 
Cf and C^'' agree with those measured for large areas of 
the 1.4 GHz BS76 survey, within the large error bars de- 
clared by Bruscoli et al. (2002); discrepancies appear in the 
comparison with results at frequencies below 800 MHz. 

• The polarization angles of the template are in good 
agreement with those measured at 2.4 GHz (D97) and cor- 
rected for Faraday rotation effects in those regions where 
the D97 position angles show a smooth dependence on co- 
ordinates. 

The last two items prove the validity of step 3. In this 
connection, we wish to stress that a perfect agreement is not 
expected at all for angular power spectra, even in the case 
of the 1.4 GHz BS76 survey. Since a conspicuous fiattenig of 
polarization power spectra is attributed to Faraday rotation, 
we expect the power spectra derived from our template to 
be somewhat steeper than those of Bruscoli et al (2002). 
This effect is not so clear due to the large error bars and 
different sky coverages, but perhaps it is already marginally 
significant. The results in Table 2 can be compared to the 
weighted averages of the 1.4 GHz angular slopes provided 
by Bruscoli et al. (2002), namely as = 1.8 ± 0.3 and as = 
1.4 ± 0.3. We note also that we are unable to find differences 
between ap and api in our template. It is an open question, 
whether such slopes will be eventually found to be equal in 
synchrotron spectra for vanishing Faraday effects. 

In conclusion, our method results in a template of the 
polarized Galactic synchrotron emission at 1.4 GHz free 
from Faraday rotation effects, which can thus be directly ex- 
trapolated to the cosmological window frequencies. In this 
range the Galactic synchrotron emission is expected to play 
the leading role in the foreground contamination of CMBP 
data. Following Platania et al. (1998), the extrapolation can 
be performed using a power law with spectral index ,3 ^ 3.0. 
As mentioned in Section 1 this template, together with the 
simulated CMBP data, provides a more reliable source map 
to test data processing and foreground separation algorithms 
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for CMBP experiments in the 20-100 GHz range. The model 
has been developed so far for a FWHM — 7° angular resolu- 
tion, matching the needs of large scale CMBP experiments 
like SPOrt, and allows to build angular spectra only up to 
£ ~ 20. At present, the position angle data represent the 
major constraint, the Heiles data being sampled on a few 
degree distance. However, we believe that the method can 
be applied at subdegree scales as well, when a complete set 
of data on these scales will be available. 
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APPENDIX A: DISTANCE L{L,B) BETWEEN 
THE SUN AND A GALACTIC HALO POINT 

The point P on the Galactic halo (radius R — 15 kpc) at 
Galactic coordinate {I, b) is the intersection between the 
sphere representig the halo and the line starting from the 
Sun position. In the Sun reference frame (see Figure Al) 
this intersection can be expressed as: 



{x - df +y^ + z^ 



cos(6) cos(i) 



cos(&) cos(i) 



cos(fe) sm(l) 



(Al) 



sin(&) 



where d is the Sun distance from the Galactic Centre (GC). 
The solution of the system provides the P coordinates 



X = dcos2(6)cos^(0 1±^/1 + 
y — xtan{l) 

z — a; tan(6)/ cos(Z) 
Therefore, the distance L(l,h) is given by 



cos-^(&) cos2(i) 



J (A2) 



L{l,b) — d cos(fe) cos(i) 



1 + 



{R^/d^ 



s'^ib) cos2(0 



(A3) 
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